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ABSTRACT 

Whilst the star formation rate (SFR) of molecular clouds and galaxies is key in understanding galaxy 
evolution, the physical processes which determine the SFR remain unclear. This uncertainty about 
the underlying physics has resulted in various different star formation laws, all having substantial 
intrinsic scatter. Extending upon previous works that define the column density of star formation 
(Esfr) by the gas column density (Egas), we develop a new universal star formation (SF) law based 
on the multi-freefall prescription of gas. This new SF law relies predominantly on the probability 
density function (PDF) and on the sonic Mach number of the turbulence in the star-forming clouds. 

By doing so we derive a relation where the star formation rate (SFR) correlates with the molecular 
gas mass per multi-freefall time, whereas previous models had used the average, single-freefall time. 

We define a new quantity called maximum (multi-freefall) gas consumption rate (MGCR) and show 
that the actual SFR is only about 0.4% of this maximum possible SFR, confirming the observed 
low efficiency of star formation. We show that placing observations in this new framework (Esfr 
vs. MGGR) yields a significantly improved correlation with 3-4 times reduced scatter compared to 
previous SF laws and a goodness-of-fit parameter = 0.97. By inverting our new relationship, we 
provide sonic Mach number predictions for kpc-scale observations of Local Group galaxies as well 
as unresolved observations of local and high-redshift disk and starburst galaxies that do not have 
independent, reliable estimates for the turbulent cloud Mach number. 

Subject headings: galaxies: high-redshift — galaxies: ISM — galaxies: starburst — ISM: clouds — 
stars: formation — turbulence 


1. INTRODUCTION 

The evolution of galaxies is driven by the intri¬ 
cate power play between interstellar turbulence and 
gravity controlling the fo r matio n of stars in giant 
mole c ular clouds (iFerrie^ 120011: iMac Low fc KlessenI 

2004 lElmegreen fc Scaldl2004 iMcKee fc Ostrikeill2007t 

Padoan et al.l I2014f) . The rate at which stars form is 
hence a pivotal quantity in tracing a galaxy’s funda¬ 
mental properties and distribution of matter and ac¬ 
tivity. Nonetheless, the functional dependence of the 
column density of star formation (Esfr) has been a 
highly debated topic for more than two decades, with 
historical parameterizations including the mean col¬ 
umn dens i ty of available gas (Eg a A (iSchmidt 19M; 
Kennicntil Il998t iBigiel et al.l 12004 iLerov et al. 200S : 

Daddi et al.l 120101: iSchruba et al.l 120111 : iRenaud et al.l 

20121: iKennicutt fc EvansI l2012h . as well as the ra- 
tio between Ep-as and the rnean freefall time i ff 
(|Krnmholz et al.Il2012l : lFederrathll2013l : lKrumhoi3l2014l l. 

However, significant scatter remains in both these ap¬ 
proaches, such that Esfr can vary by more than an order 
of magnitude for any given input Egas or Egas/tff • Follow¬ 
up analyses of these theoretical models via computer sim¬ 
ulations have determined that the observed scatter may 
be primarily attributed to the physical variations in the 
sonic Mach nu mber (A4) of th e turbulence in the star¬ 
forming clouds (jFederratbll2ni3[l . This results in separate 
Al-dependent relations for the SFR, which cover the ob¬ 
served range of Esfr- The aim of this work is hence to 
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unify these separate laws and develop a unique, universal 
relation. 

In Section [2] we introduce and derive our new SF law 
based on the PDF and Mach number of interstellar tur¬ 
bulence. Section [3] describes the observational sample 
utilized to assess our new SF law. In SectionUwe present 
our findings and compare our new model to previous pa¬ 
rameterizations. In Section [5] we invert our new SF law 
to make sonic Mach number predictions for extragalactic 
sources. Finally, Section [6] summarizes our conclusions. 

2. A NEW STAR FORMATION RATE DESCRIPTOR 

Here we derive an improved SF law, which takes into 
account the distribution of gas densities (PDF) present in 
galactic clouds and the interstellar medium (ISM). This 
follows from the finding that the observed scatter in the 
SF law can be principally attributed to the variation in 
the son ic Mach numbe r of the clouds and galaxies ob¬ 
served (|Federratbll20I3ll . 

2.1. The probability density function (PDF) of 
interstellar gas 

In 2012 Krumholz, Dekel and McKee initiated the re¬ 
finement of the traditional Kennicutt-Schmidt law by 
dividing the gas surface density by the freefall time of 
the gas (Egas/tff), indeed achieving a stark improve¬ 
ment in star formation rate correlation (jKrumholz et 
I20I2I: iFederrathI I20I3I: iKrnmhol j I20I4|1 . However, be¬ 
cause this model implements only the division between 
the average surface density and the average freefall time 
of the gas, important information regarding the wide 
distributio n of densities within the ISM and molecu- 
lar clouds (iKainulainen et all2009l: [Schneider et al.ll201^ 
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l2013t iKainulainen et aT1l2013L l2014fl is not taken into ac¬ 
count. We hence extend upon the important previous 
work by iKrumholz et al.l (l2012f) to incorporate the den¬ 
sity PDF, with the ultimate goal of eliminating the re¬ 
maining scatter in the previous relation. 

In our new SF law, we must consider the 
fact that denser gas forms stars at a highe r rate 
(IKrumholz fc McKed 120051 : iPadoan fc Nordlundl 120111 : 
lHopkinjl2013[l because gas with higher density p has a 
shorter freefall time than more diffuse gas, 

“ \lS2Gp’ 

where G is the gravitational constant. Note that 
the concept of a mult i -freefall tim e wa s pioneered by 
iH ennebelle fc Chabrieil (|201lL 1201311 and iChabrier et al.l 
()2014ll . It suggests that, due to the clumpy nature 
of molecular clouds, the typical timescale for star for¬ 
mation is not the average timescale of the clouds, but 
the density-dependent timescale of each collapsing sub¬ 
structure within the clouds. This multi-freefall concept 
h as recently received support f rom numerical simulations 
iniFederrath fc KlessenI (1201211 and from observations in 
lEvans et al.1 (1201411. _ 

In contrast. IKrumholz et al.l (j2012l l and iFederratdil 
(|2013ll had only correlated Esfr with Egas(po)/te(Po) 
(hereby denoted as (Egas/Osingie-ff); he. only at a sin¬ 
gle freefall time evaluated at the mean density po- It 
is conceivable, however, that the true correlation must 
be in the form of an integral over the whole density 
PDF, J'[Egas(p)/tff(p)]ppdp. In the following, we com¬ 
pute this integral by making the standard assumption 
that the PDF follows a log-normal distribution as the 
initial condition for star formation, 

where the density is written in logarithmic form, normal¬ 
ized to the mean density 

s = \n{p/po). (3) 

This transformation of variables from p to s has ad¬ 
vantages during the integration step below and allows 
us to write the PDF, Equation ([2), in its standard log¬ 
normal form. The mean logarithmic density sn is related 
to th e variance by sn = —o'il2 (iVazauez-Seniadenil 
199411. As d e rived by iPadoan fc Nordlundl (|2011ll and 
Molina et al.l ()2ni2f l. the logarithmic density variance is 
given by 

a2=ln(^l + 6W^), (4) 


parameterized by the sonic Mach number AI, the tur - 
bulent driving parameter h (|Federrath et al.l[2008Ll2010ll . 
and the ratio of thermal to magnetic pressure, plasma /3. 

The log-normal PDF, Equation ([2]), and the associated 
density variance - Mach number relation. Equation (|^ . 
are observed to provide a good approximation to the real 
PDF i n molecular clouds (IKainulainen e t al.ll2009l: iBrunt 
201C : iSchneider et al.l 120121 120131: Kainulainen et al. 


2014 1. in the Galactic c entre ([Rathborne et al. | 


on large Galactic scales ()Berkhuiisen fc Fletche: 


2014 1. 


2001 , 


and even in extra-galactic sy stems (iHughes et al.l[2M^ 
iBerkhuiisen fc Fletcher! [2 01 31 . 


2.2. Derivation of the maximum (multi-freefall) gas 
consumption rate (MGCR) 

For a rigorous derivation, we start by defining a carte¬ 
sian coordinate system in which the line-of-sight toward 
the cloud or galaxy is in the z-direction and maps of col¬ 
umn density and star formation rate surface density are 
in the xy-plane. First, the gas column density is defined 

as Egas(p) = /q pdz, where is the scale height of 
the cloud or galaxy for which Egas is to be determined. 
From observations and simulations discussed in the pre¬ 
vious section, we know that Egas and p follow closely 
a log-normal distribution. Equation ([5]). Thus, in order 
to compute an average of a density-dependent variable, 
we simply have to integrate that variable over the en¬ 
tire PDF. For the average gas column density, for exam¬ 
ple, we would evaluate f EgasPsgasdSgas = // pppdpdz = 
ff po exp(s)p(s)dsdz = po f dz = poH^, which is in¬ 
deed the average column density Egas(po)- Now we 
follow exactly the same mathematical procedure, but 
for the combined density-dependent variable Egas/tff = 
Sgas (Po)/iff(po)exp(3s/2), i.e.. 


(5Igas/i) 


multi—ff 


/O 

Scr, 


^gas(p) 

ts{p) 

.ifo) r 


te(po) 


Ppdp 

0 

{plPof^^Ppdp 


^(Po) 


iff(po) J-c 


exp 


-s ) p{s)ds 


^gas 


(po) 


tsipo) 


■ exp -cr; 


(5) 

( 6 ) 

(7) 

( 8 ) 


Note that in the second step, we use the exact scaling 
Sgas/iff P^^^, because Egas ~ P and te The 

third step transforms variables from p to s via Equa¬ 
tion 0 and via the identity Ppdp = p{s)ds. The last 
step. Equation ([8|), is the analytic solution of the inte¬ 
gral over all densities. 

We call this new quantity derived above, the maximum 
gas eonsumption rate or multi-freefall gas consumption 
rate (MGGR), denoted (Egas/t)muiti-ff- 


2.3. The single-freefall to multi-freefall correction factor 

We see that Equation ([8|) includes the multi-freefall 
correction factor given by exp(3(Ts/8), which corrects the 
single-f reefall average Pa] /iff ( pn) = (Ep-as/t)sinH e-ff 

used in IKrumholz et akl (120121 1 and iFederrathI (120131 1 for 
the underlying PDF of densities. 

We can now insert the density variance from Equa¬ 
tion (|T]) into Equation (|5]) and obtain the multi-freefall 
correction factor. 


Gmulti—ff 


(Egas/t)niulti—ff 
(^gas/t)single—ff 


(^1 + b'^M^ 



3/8 


(9) 

Figure ^ shows the (Egas/t)multi—ff/(Sgas/0single—ff 
correction factor. Equation as a function of the 
Mach number A4 for three values of the turbulent driv¬ 
ing parameter b and—for simplicity—neglecting mag- 
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Table 1 

Summary of observational sources considered in deriving and testing our new SF law 


Data Source 

Source Type/Name 

logic ^SFR 

loglO ^gas 

logio(^) 

single—ff 

M 

Range 

Correction 

Factor 

logio(^) 

multi—fF 

Bolatto et al. ('2011') 

Small Magellanic Cloud (SMC) 

-2.15 

1.10 

-0.670 

16-200 

12.6 

0.429 

Yusef-Zadeh et al. f2009') 

Central Molecular Zone (CMZ) 

-0.750 

2.08 

0.900 

50 

9.47 

1.88 

Heiderman et al. ('2010') 

Taurus 

-0.750 

2.04 

1.60 

10 

1.11* 

1.64 

(HIO) 

C2D -1- GB Clouds 

-0.135 

1.89 

1.68 

5-20 

2.89 

2.14 


Class I YSOs 

0.208 

2.24 

2.40 

1-2 

1.06 

2.43 


Class I YSOs (upper limit) 

-0.147 

2.01 

2.15 

1-2 

1.06 

2.17 


Flat SED YSOs 

0.332 

2.28 

2.47 

1-2 

1.06 

2.50 


Flat SED YSOs(upper limit) 

-0.090 

2.00 

2.10 

1-2 

1.06 

2.13 

Wu et al. (20101 (WIO) 

HCN(l-O) Clumps 

1.17 

2.95 

3.46 

2-5 

1.40 

3.60 

Gutermuth et al. (2011^ 

Class II YSO counts in 







(Gll) 

eight molecular clouds 

-0.787 

1.64 

1.04 

5-20 

2.89 

1.50 

Lada et al. ('2010^ 

Molecular clouds at Aj^ >0.1 

-0.941 

1.46 

0.722 

5-20 

2.89 

1.18 

(LIO) 

Molecular clouds at A]^ > 0.8 

0.545 

2.37 

2.57 

5-20 

2.89 

3.03 


Notes. Column 1: Reference. Column 2: Source classification. Column 3: Geometric mean of Ssfr- Column 4: Geometric mean of Sgas- 
Column 5: Geometric mean of (Sgas/t)single—ff • Column 6: A4 ranges whose geometric mean of the upper and lower limit were used to 
calculate the average Mach number. Column 7: Correction factor based on the average Mach number as attained in Equation Note 
that we assume plasma /3 —>■ oo (i.e., no magnetic field c orrect ion) for all except Taurus’ correction factor (see ★), for which we use the 
available measurement of /3 = 0.02 from fHever Brunt! 1)20121 . Column 8: Our new SFR descriptor, (Sgas/t)multi—ff ? called maximum 
(multi-freefall) gas consumption rate (MGCR), computed via Equation and derived in Equation I®- 



Figure 1. The multi-freefall correction factor derived in Equa¬ 
tion lO as a function of the sonic Mach number for three values of 
the turbulent driving parameter b whilst neglecting magnetic fields 
(/3 —>• oo). The values in the figure represent the correction factor 
are applied to existing measurements of (Sgas/i)single—ff order 
to correct them to (Sgas/i)multi-ff • 


netic fields, i.e., /3 — >■ oo. We see that the multi- 
freefall correction is always positive and increases with 
increasing Mach number. The correction factor is about 
an order of magnitude for typical cloud Mach num¬ 
bers, A4 ^ 5-50, which means that our correction to 
the previous SF laws will be significant. In the follow¬ 
ing, we will assume a fixed turbulent driving parameter, 
b = 0.4, representing a natural mixture. However, we 
emphasize that b is not likely fixed ac r oss all molecu¬ 
lar clouds (lBruntll2010HPrice et al.ll2011l : iGinsburg et al.l 
120131 : iKainulainen et al.ll2014D . but in the absence of a 
direct measurement of b, we presume a standard value 
for all clouds. 


3. OBSERVATIONAL SAMPLE SELECTION 

We selected our sample of 11 galactic objects and 
1 extragalactic object based on the availability of gas 
or dust column density (Sgas) and SFR column den¬ 
sity (Ssfr)- Submillimeter observations are the pri¬ 
mary source for the study of molecular gas clouds 
in galaxies. In combination with young stellar ob¬ 
ject (YSO) counts, infrared and ultraviolet luminosi¬ 
ties, Ssfr measurements may be obtained. As our 
new SFR descriptor is primarily dependent on the 
sonic Mach number of the gas in each observed source, 
availability of such estimates were vital in selecting 
relevant dat a. Mach number est i mates for Milky 
Way clouds (iHeiderman et al.l 120101 iLada et aP 120101: 
Gutermuth et al.l l2011h. molecular clumps (iWu et al.l 
201011 . YSO ([Heiderman e t al.l l20 inli and the Central 
Molecular Zone ICMZl (jYusef-Zadeh et al.Jl2009ll as well 
as the Small Magellanic Cloud ISMCi (iBolatto et al. 

m were taken from and are summarized in iFederratb 
. A summary of the data and sonic Mach number 
estimates utilized are listed in Table [TJ 
Derivation of (Egas/t)muiti-ff was calculated via Equa¬ 
tion ([9]) using b = 0.4 and /3 —>■ oo for all cases except the 
Taurus molecular cloud. Ta urus is known to exhib it high 
levels of magnetic activity (iHever fc Bruntl[2012h . so an 
estimate of /3 = 2A4\/A4^ ~ 0.02 (with the sonic and 
Alfven Ma ch number, M = 10 and A4^ = 1, respec¬ 
tively; see iHever fc BruntI 1201211 was applied to obtain 
Taurus’ multi-freefall correction factor. 

Previous SF laws have reported data points at¬ 
tained from individual clouds and assigned a linear 


correlation between Sc-fr and Sp-as or (S 


Jt 


gas / ^ single—ff 

(jKennicutt fc Evansll2012t iKriimholz et al1l2012l l. How¬ 
ever, because we do not have sonic Mach number mea¬ 
surements for individual clouds, utilizing the data points 
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Table 2 

Summary of the results of the statistical tests applied to the data 
and the fits generated to model their relationship 


logjQ SsFH, correlator 

Offset 

Slope 

Scatter 



Two-parameter fits 



^•^§10 ^gas 

-3.70 

1.70 ± 0.21 

3.83 

0.864 

(^gas /^)single —ff 

-1.61 

0.78 ±0.05 

1.07 

0.962 

^^§10 multi —ff 

-2.37 

1.01 ±0.06 

0.99 

0.965 


One-parameter fits 



(^gas /^)single —ff 

- 1.99 

1 (fixed) 

3.34 

0.881 

^^§lo(^gas A) multi —ff 

-2.34 

1 (fixed) 

1.00 

0.965 


Notes. Column 1: SFR descriptor in log space. Column 2: Offset 
of the relation’s robust line fit in log space. Column 3: Slope of the 
relation’s robust line fit in log space. Column 4: Relative scatter 
of the data. Column 5: Goodness-of-fit, value = 1 would 
indicate a perfect fit). 

from each cloud to derive a law from our new descrip¬ 
tor would not be an accnrate method of derivation. We 
can only obtain approximate values of sonic Mach num¬ 
ber for populations of clouds, so in order to arrive at a 
fair comparison for each relation we report the geometric 
mean values of Ssfr, Sgas, (Sgas/i)singie-ff, and M for 
each group of clouds and apply robust line fits to these 
data points. The error bars associated with each point 
were derived from calculating the standard deviation of 
the mean. Since the previous SFR descriptors Egas and 
(Egas/t)singie-ff do not depend on the sonic Mach num¬ 
ber of the clouds, uncertainty in the values of A4 was 
not propagated when deriving values of (Egas/t)muiti-ff 
in order to objectively evaluate differences between our 
new SFR descriptor and previous works. 

4. RESULTS: AN IMPROVED STAR FORMATION LAW 

Figure[2]shows a direct comparison of previous SF laws, 
Esfr vs. Egas (left panel) and Esfr vs. (Egas/t)singie-ff 
(middle panel), together with our new MGCR corre¬ 
lator, Esfr vs. (Egas/t)muiti-ff (right panel), includ¬ 
ing all the data listed in Table [T] For the classical 
Kennicutt-Schmidt relation (left-hand panel), we apply 
a two-parameter robust power-law fit, with the offset 
and the slope of the power law as free parameters. For 
the iKrumholz et al.l (1201211 relation (middle panel) and 
for our new multi-freefall relation (right-hand panel), we 
respectively apply a one-parameter robust fit with the 
offset as the only free p arameter (the slope is fixed to 
unity). Thus, both the IKrumholz et al.l (1201211 single- 
freefall SF law and our new multi-freefall SF law predict 
a direct linear relationship between the actual SFR and 
the maximum possible single-freefall or maximum possi¬ 
ble multi-freefall gas consumption rate, respectively. We 
also apply the two-parameter robust fit to both these re¬ 
lations in order to gauge how much they deviate from 
the assumed unity slope. 

In order to quantitatively evaluate the relative good¬ 
ness of fit for each relation, we have applied two sta¬ 
tistical tests to the data generated by each parameter¬ 
ization and their corresponding lines of best fits, sum¬ 
marized in Table [21 The first value, which we define as 
the relative scatter in each relation, is the unweighted 
value normalized with respect to our new SFR descriptor, 
(Egas/t)muiti-ff ■ Scatter values greater than unity hence 


indicate a greater degree of scatter than that present in 
our new relation. The second test applied attains a value 
for the goodness-of-fit, The value is the coeffi¬ 
cient of determination and has a range of 0 to 1. An 
R^ = 1 indicates a model that perfectly fits the data. 

Both the statistical tests we have applied show 
a significantly tighter correlation between Esfr and 
( Egas/t) multi—ff than with either (Egas/t)singie—ff or Egas, 
as summarized in Table |2j For the (Egas/t)singie-ff and 
(Egas/t)muiti-ff parameterizations, this conclusion is fur¬ 
ther emphasized by the behavior of their residuals and 
gradients of their two parameter least-squares fits. The 
two-parameter fits suggest that there is significant in¬ 
trinsic correlation remaining in the (Egas/t)singie-ff pa¬ 
rameterization, for which we obtain a best-fit slope of 
0.78 ±0.05 instead of unity. In contrast, our new MGCR 
correlator yields a slope of 1.01 ± 0.06, consistent with 
an intrinsic slope of unity, i.e., a truly linear correlation 
between Esfr and (Eps/t)muiti-ff- 

Our results in Figure [2] and Table [2] strongly 
suggest that the remaining s catter in the p r evious 
Esfr vs. (Egas/t)singie-ff law by IKrumholz et al.l (l2012f) 
can indeed be primarily attri buted to systema tic vari¬ 
ations in A4 as suggested by iFederratiil (j201,‘lfR . The 
accounting of said Mach number variations in our new 
SF law has eradicated the discrepancies of up to an or- 
der of magnitud e , colla psin g the scatte r ed rel ations from 
IKrumholz et al.l (j2ni2| l and iFederratbl (1201,811 into a sin¬ 
gle, more universal SF law, 

F^SFR — 0.4% X (Fjgas/^)multi— ff 

= 0.4% X (Sgas/^)single—ff ^ C^multi—ff 

/ R 

= 0.4% X (EgasA)single-ff X fl ± 

( 10 ) 

The reduction in scatter in this SF law strongly indicates 
that our new model is a physically meaningful descriptor 
of the SFR. 

Thus, our new SF law given by Equation m pre¬ 
dicts that the SFR is equal to ~0.4% of the MGCR in a 
molecular cloud or galaxy. This value is slightly smaller 
than the 1% value of th e mean gas cons u mptio n rate, 
(Egas/t)singie-ff, uscd in iKrumholz et al.l (1201211 . The 
reason for this is that the MGCR is naturally larger 
than the mean gas consumption rate, because we take 
the density-dependent freefall time into account in our 
model. 

5. APPLICATION OF THE NEW MODEL: MACH NUMBER 
PREDICTIONS FOR EXTRAGALACTIC SOURCES 


The IKrumholz et al.l (|2012ll SF law had util ized extra- 


galactic data of disk and starburst galaxies (iKennicutt 

1998 

Bouche et aLll2007l: 

Daddi et al.ll2010l:IGenzel et al. 

201 c 

Tacconi et al.l 120 1C 

) in addition to Milky Way ob- 


servations. However, it is difficult to accurately gauge the 
values of A4 for extragalactic sources as we are not able 
to well resolve individual molecular clouds, especially at 


^ S ee als o IKrumholz fc McKed Il200.‘]ll. IPadoan fc Nordlundl 
([201JJL__ [HennebeUg^^^CTiabrTCT 1 I 2 OI 1 I [20T^ and 
IFederrath &: Klessenni20T^ on the detailed derivation of how the 
SFR depends on the Mach number in the single-freefall versus the 
multi-freefall formulation. 
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Figure 2. Bottom panels: Ssfr versus Sgas (classical Kennicutt-Schmidt relation; left), Ssfr versus (Sgas/f(single—ff UKrumholz et aO 
120121 . middle), and our new model, SgFR versus (Sgas/f)multi—fl (right) for the observational data in Table [T] Power-law fits are shown as 
solid lines and the respective goodness-of-fit parameter (where = 1 corresponds to a perfect fit) is shown in each panel, as well as 
the normalized scatter. We also apply a two-parameter robust line fit to these relations, shown as the dotted lines. We find that our new 
multi-freefall SF law provides the best correlation with a reduced scatter by a factor of 3-4 compared to the previous SF laws. Top panels: 
The residuals of each relation in log space versus each of the parameterizations being analyzed. The solid and dotted lines in each panel 
are respectively the fits shown in the bottom panels. 


Table 3 

Mach number predictions for extragalactic systems obt ained by 
inverting our improved SF law, Equation lIlQ]l 


Data Source 

Redshift 

Galaxy Type 

A4 Estimate 

Kennicutt fl998) 

0 

Disk 

4.0ti° 


~ 0 

Starburst 

10 + 21 

Daddi et al. ('2010') 

1-3 

Disk 

letu 

Tacconi et al. ('2010') 

1-3 

Disk 


Genzel et al. ('2010') 

1-3 

Disk 

K 7 -I-O .5 


1-3 

Starburst 

Cl+120 

^-*^-36 

Bouche et al. ('2007') 

1-3 

Starburst 

71 +120 

' ^-45 


Notes. Column 1: Reference. Column 2: Redshift range. Col¬ 
umn 3: Galaxy type. Column 4: Geometric mean A4 prediction 
from inverting the improved SF law. Equation m, as well as the 
upper and lower limits of this A4 estimate. 

high redshifts. A4 estimates require measurements of 
the sound speed and turbulent velocity dispersion of a 
cloud, which in turn requires high-resolution molecular 
line data and a temperature measurement. Such data 


are hard to obtain for molecular clouds in external galax¬ 
ies, as a complete census of the star-forming molecular 
cloud population in a given galaxy is needed at suffi¬ 
ciently high resolution to capture the cloud-scale Mach 
number. There are only very few studies of nearby galax¬ 
ies (such as M51) w hich start to resolve molecular clouds 
in external galaxies (jHughes et al.l[2r)T^ , but even those 
do not necessarily yield a complete census of clouds and 
it is not clear that all cloud properties are converged at 
the telescope resolutions available to date. 

We were hence unable to include extragalactic sources 
that lack A4 estimates in the deduction of our new SF 
law. They were similarly omitted in our analysis of pre¬ 
vious laws in order to arrive at a fair comparison be¬ 
tween the different parameterizations of Ssfr in Figure[2] 
and Table [21 Instead, we invert our new model, Equa¬ 
tion m, in order to compute A4 estimates for these ex¬ 
tragalactic sources; predictions to be tested by future ob¬ 
servations of resolved CO maps in extragalactic sources. 
Our A4 predictions are summarized in Table [31 Since the 
presence of vigorous star formation in starburst galaxies 
would result in more violent turbulence of gas compared 
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to disk galaxies, it is expected that the starburst galax¬ 
ies have greater A4 than their disk galaxy counterparts of 
the same redshift range, which is indeed the systematic 
trend we find in Table |3l This trend is consistent with 
the analyses and conclusions presented in IChabrier et alJ 

dS). 

6. CONCLUSIONS 

We derived a new SF law given by Equation OT, which 
is based primarily on the sonic Mach number of turbu¬ 
lence, as a result of the density PDF of molecular, star¬ 
forming gas clouds. We find that the SFR is equal to 
^ 0.4% of the maximum (multi-freefall) gas consump¬ 
tion rate (MGCR), which we derived in Section [2J We 
compared our new model to previous parameterizations 
of the SFR and determined quantitatively that our new 
SF law provides a tight linear relation between Esfr and 
MGCR, with a factor of 3-4 less scatter compared to any 
previous SF law. By inverting Equation (nni), we predict 
the Mach numbers (At) of the star-forming molecular 
clouds in extragalactic sources. Our predictions are sum¬ 
marized in Table |3l which anticipate testing via future 
submillimetre observations. 
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work. C.F. acknowledges funding provided by the Aus¬ 
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fully acknowledges the support of an ARC Future Fel¬ 
lowship and ARC Discovery Project DP130103925. 
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